! { dg-do compile }

      SUBROUTINE PROPAGATE(ICI1,ICI2,I,J,J1,ELEM,NHSO,HSO
     *                    ,MULST,IROOTS)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      COMPLEX*16 HSO,ELEM
      DIMENSION HSO(NHSO,NHSO),MULST(*),IROOTS(*)
      ISHIFT=MULST(ICI1)*(I-1)+1
      JSHIFT=MULST(ICI2)*(J-1)+1
      DO 200 ICI=1,ICI1-1
         ISHIFT=ISHIFT+MULST(ICI)*IROOTS(ICI)
  200 CONTINUE
      DO 220 ICI=1,ICI2-1
         JSHIFT=JSHIFT+MULST(ICI)*IROOTS(ICI)
  220 CONTINUE
         DO 150 MSS=MS,-MS,-2
           IND1=ISHIFT+K
           IND2=JSHIFT+K
           HSO(IND1,IND2)=ELEM
           HSO(IND2,IND1)=DCONJG(ELEM)
  150    CONTINUE
      END


